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Abstract 

We generalize the collision term in the one-dimensional Boltzmann-Nordheim 
transport equation for quasiparticles that obey the Haldane exclusion statis- 
tics. For the equilibrium situation, this leads to the "golden rule" factor 



^ ' for quantum transitions. As an application of this, we calculate the density 



response function of a one-dimensional electron gas in a periodic potential, 
assuming that the particle-hole excitations are quasiparticles obeying the new 
statistics. We also calculate the relaxation time of a nuclear spin in a metal 
using the modified golden rule. 
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1. INTRODUCTION 

Recently, Haldane has proposed a new exclusion statistics for confined quasiparticles, 
that interpolates between the Fermi and Bose quantum behaviours, and gives rise to partial 
blocking of a singly occupied state. *^^) There has been considerable theoretical activity con- 
cerning study of statistical and thermodynamical properties of systems of particles obeying 
the generalised exclusion statistics. Examples of quasiparticles obeying this new statistics 
include spinous in antiferromagnetic spin chain with inverse square exchange interaction'-^-*, 
and two-dimensional anyons restricted to the lowest Landau level in a strong magnetic 
field. ''^-'~''^-' Other examples are the Calogero-Sutherland model*^^^"*^^^ in one dimension, and 
particles interacting with a two-body delta-function potential in two- dimensions in the ther- 
modynamic limit. *-^-' Haldane *^^^ has analysed a Heisenberg spin chain in one-dimension with 
inverse- square coupling, and shown that under certain conditions the quasi-particle excita- 
tions obey a particularly simple form of the generalised statistics. Despite all these theo- 
retical studies, it is not entirely clear under what conditions a physical or chemical system 
will experimentally manifest this generalised statistics. The best evidence so far comes from 
recent neutron inelastic scattering experiment ''^°-' on the compound KCuF^ which is a one- 
dimensional Heisenberg antiferromagnet above 40^K. The observed inelastic scattering is 
best fitted by spinon excitations in a spin chain whose interactions fall off as the inverse 
square of the lattice distance*^^^^ . The dynamic spin correlation function for such a system 
has been calculated by Haldane and Zirnbauer.'^^^) 

The present paper addresses a question of current interest, namely how to calculate 
transport properties of particles obeying the generalised exclusion statistics. Specifically, 
how does the Boltzmann transport equation get modified in such a situation? It may be 
recalled that the collision term in the Boltzmann equation was modified by Nordheim'-^^'' 
to incorporate the effect of quantum statistics. Account was taken of the Pauli blocking 
of the final state in binary collision events for particles obeying Fermi statistics. On the 
other hand, an enhancement factor was introduced (as in the Einstein coefficient of induced 



emission'^^^) ) in the final state occupancy for Bose statistics. Tlie suppression and enliance- 
ment factors, (1 — /) and (1 + /), appropriate for Fermions and Bosons, respectively, are 
routinely used in many-body calculations. What are the corresponding factors for the sys- 
tems of particles obeying the new statistics? We have presented the answer to this question. 
Future calculations of transport properties of such systems will find these results relevant. 
Besides generalizing the Boltzmann-Nordheim equation, we have also generalized the Fermi 
Golden Rule. 

The plan of this paper is as follows. In section 2, we examine changes in the structure of 
the collision term in the Boltzmann equation due to the incorporation of the new statistics. 
This yields the "golden rule" factor for quantum transitions. As an application of this, 
we examine, in section 3, the response of a one-dimensional electron gas to an external 
periodic potential, and calculate its density response function. Whereas the ground state of 
the electron gas is taken to be the conventional Fermi system, we assume that the particle- 
hole excitations are quasiparticles obeying the new statistics. The magnitude of the density 
response function is then shown to be appreciably enhanced, thereby lowering the Peierls 
transition temperature. In section 4, we generalize the Fermi Golden Rule for quasiparticles 
obeying the new statistics and apply it to the calculation of a relaxation time of a nuclear 
spin in a metal. Both the applications are schematic, with a view to explore qualitative 
changes from the standard results. 



2. GENERALIZED BOLTZMANN-NORDHEIM EQUATION 

We now derive the Boltzmann equation (specifically, the collision term) for particles 
obeying Haldane's statistics. For the derivation of the collision term, we make the usual 
assumptions of (i) the dominance of binary collisions, (ii) the hypothesis of molecular chaos 
("stosszahlansatz"), and (iii) the slow variation of the distribution function f{x,p,t) over 
distances and times of the order of characteristic interaction lengths and durations, respec- 
tively. Recall, f{x,p,t)dx dp/{2Ti) is the number of particles in the phase-space volume 
element dx dp at time t. (We set /i = 1 throughout this paper.) We wish to derive an 
expression for the collision term C in the Boltzmann equation: 

dt dx dp 

C is the rate of change of / caused by collisions. Collisions of the type p'p'i —^ ppi tend to 
increase the population of particles having momentum p, while those of the type ppi -^ p'p\ 
tend to decrease it. Let w{p'p'i -^ ppi)dp dpi be the transition probability per unit volume 
and per unit time that two particles having incoming momenta p' and p[ are scattered 
with outgoing momenta in the ranges {p,p + dp) and {pi,Pi + dpi). Then according to the 
hypothesis of molecular chaos, for classical (Maxwell-Boltzmann (MB)) particles the net 
increase in the number of particles in dx dp due to collisions occurring during dt is*^^^-* 

/ f'dp.f[dp[.w{p'p[ -^ ppi)dp dpi dx dt 
—f dp.fidpi.w{ppi — > p'p\)dp' dp\ dx dt 

where /' = f{x,p',t), etc., and the integration is only over dp'dp[dpi. Using the detailed- 
balance property w{p'p[ -^ ppi) = w{ppi -^ p'p'i) and inserting a factor | to take into 
account the fact that a final state with momenta (p',Pi) is indistinguishable from that with 
momenta {p'i,p') we get 

C = \jdp,dp'dp[ w{pp, ^ p'p[){f'f[ - fh). (2) 



For particles obeying quantum statistics, however, the above derivation needs to be mod- 
ified because now the transition probabihty depends also on the occupancy of the final state. 
Thus, if'fi - fh) in (2) should be replaced by /'/{ F{f)F{f,) - f f\ F{nF{f[), where F 
takes account of the effect of the quantum statistics on the accessibility of the final state. 
For example, in the case of Fermi-Dirac (FD) statistics, F{f) = (1 — /) to incorporate Pauli 
blocking, and for Bose-Einstein (BE) statistics, F{f) = (1 + /) to take care of Einstein 
enhancement. According to the ilf -theorem in statistical mechanics, the entropy production 
vanishes if and only if 

rnF{f)F{f,) = ff,F{f)F{f[). (3) 

Solution of this equation provides the shape of the equihbrium distribution function f^q. 
Conversely, as we proceed to show, if f^g is known, (3) can be used to deduce the functional 
form F{f). Rewriting (3) in the form 

/ /i f f[ 

F{fyF{h) FifYFifiY 

taking logarithms of both sides, and using the fact that energy (E) is the appropriate 
summational invariant in the collision, we get 

ln[//F(/)] = -a(E-6), 

where a and b are constants to be identified with inverse temperature (T) and chemical 
potential (/i), respectively. Thus, we get finally 

F{f) = exp[a{E - b)]f = ({E)/, ({E) ^ exp[(E - fi)/T], 

where we have set the Boltzmann constant to unity. For the MB, FD and BE statistics, 
the equilibrium distribution functions / are (~^, l/((^ + l) and l/((^ — 1), respectively. Hence, 
the corresponding F's are 1, (1 — /) and (1 + /), as expected. 

Now, for particles obeying the Haldane statistics, the equilibrium distribution function, 
as shown by Wu'-^^ is 



f{E) = l/{u{E) + a), (4) 

where a (0 < a < 1) is the statistical interpolation parameter (a = corresponds to Bosons 
and q; = 1 to Fermions) and uj{E) satisfies 

o;(E)"[l + cu(E)]i-" = exp[{E - ^)/T]. (5) 

It is straightforward to show that, in this case, 

F{f) = {l-afr[l + {l-a)fY'^. (6) 

In particular, for semions (a = |j we get 

F(/) = (l-//2)"=(l + //2)"^ (7) 

Thus, the Boltzmann equation (1) for particles obeying Haldane's statistics is given by 

dt dx dp 

= ^Jdp,dp'dp[ w{pp,^p'p[)[f'f[ F{f)F{f,)-ff, F{f')F{f[)], 

where F is given by (6). 



3. DENSITY RESPONSE FUNCTION 

In calculating the density response function xiQ) of the one-dimensional electron gas, we 

follow the treatment given by Kagoshima et alM^^ For a periodic potential V{r) = ^ Vqc**^'', 

Q 

the charge distribution of the electrons in a length L is distorted by 5p{r) = -f y^ppe^'^'', 

Q 
and x{Q) is just the linear response function, obtained from 

PQ = -Vq x{Q), 

by calculating Sp{r) perturbatively. Consider an electron in the initial state k, absorbing a 
phonon of wave number Q {—Q), reaching the final state k + Q {k — Q). Then, using FD 
statistics, the usual expression for pg is*^^^^ 



pq - ^ 2^ 



/fc(l — fk+Q) /fc(l — fk-Q) 



Ek — Ek+Q Ek — Ek-Q 



(8) 



k 

where we have neglected the phonon energy Uujq in comparison to the electron energy 
Ek = h A;^/2m, and A^ is the number of electrons in length L. The dummy variable k in the 
second term may be replaced by {k + Q), resulting in the cancellation of the bilinear terms 
in /'s, and yielding 

X(Q) = ^E#±^- (9) 

At temperature T = 0, we have fk = 1 and fk+q = 0. In one dimension, this results in a 
logarithmic singularity for transitions from k = ^.kp to k = =F^f with Q = ^2kF- This is 
a manifestation of the well-known Peierls instability. At nonzero temperature T = (3~^^ the 
singularity is replaced by a finite result, as may be verified by a straight-forward calculation: 

(n 01- \ :D(Ep) p/3/2 tanhx 
X[Q = 2kF) = ^.. / dx, 10) 

where D{Ef) is the density of states at the Fermi level and the energy integration is carried 
out only in the range \E~Ef\ < e^ ^ Ef- We repeat this calculation for x by replacing the 
factors /(I — /) in (8) by the new factors fF{f) derived in section 2. There is of course, no a 
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priori justification for assuming tliat tliese excitations obey tlie new statistics. Nevertlieless, 
if tliey did, what would be the experimental signals in a one-dimensional system? The 
enhancement of x{Q) near Q = 2kF has a direct bearing on the estimation of the Peierls 
metal-insulator transition temperature, as well as on the Kohn anomaly. It may therefore 
be a worthwhile exercise to reexamine x{Q) with the modified golden rule. 

It is simplest to illustrate the calculation with semions {a = 1/2), where the equilibrium 
distribution is explicitly known (see (7)). The expression in the square brackets in (8) is 
then replaced by 



1/2 . /. 1.2 \l/2' 



fk (l - ifk+Q) fk (l - i/j 



Ek — Ek^Q Ek — Ek^Q 



(11) 



Note that on changing the dummy variable k ^ k + Q in the second term in (11), there is 
no cancellation of the bilinear terms in / as earlier. However, at T = 0, the quasiparticle 
excitations near the Fermi surface have fk = 2, fk+q = 0, and hence, there is an overall 
enhancement of x{Q) by a factor of 2, again with a logarithmic singularity at Q = 2kF- 
At nonzero temperature, the semionic distribution function / given by (7) may be used to 
evaluate x{Q)^ and one obtains after straight-forward algebra 

X[i^-ZKF)- 2Ar ^io a; (l + 4e-4-)i/2(i+4e4x)i/2' ^"^^^ 

which should be compared with (10). The integrals in (10) and (12) are plotted as a function 
of €3/3/2 in Fig. 1, and again show an enhancement of x{Q)> but by a factor somewhat 
greater than 2. Although there is an enhancement in the response function in our simple 
model, because the temperature dependence of it does not change, it may not be possible 
to detect this effect experimentally. 



4. GENERALIZED FERMI GOLDEN RULE AND RELAXATION TIME 

We now consider generalization of the well-known Fermi Golden Rule, for particles 
obeying the fractional exclusion statistics. Consider a two-particle scattering process 
k + p ^ k' + p', in one dimension. Let the two particles be distinct. We have in mind 
processes like the electron-nuclear interactions in metals. According to the Fermi Golden 
Rule, the transition probability per unit time for such a process is given by 

Wi^f = 27i\ < i\Hint\f > \^S{Ei - Ef), 

where Hint is the interaction Hamiltonian and the delta function ensures energy conservation 
between the initial and the final states. Typically in electron-nuclear interactions (as in 
magnetic relaxation in solids) one can neglect the nuclear recoil to a very good approximation 
and therefore one can replace S{Ei — Ef) by 6{Ek — E^r), k and k' being the initial and final 
momenta of the scattered electron. We work in this limit to keep the calculation simple and 
also because this is the relevant limit for many physical applications. The total transition 
probability per unit time at nonzero temperature is then given by*^^^^ 

W^^f = 2njdkdk'\ < t\H,nt\f > \^S{Ek - Ek')p{k)p{k')f{Ek)F{Ek'), 

where p{k) denotes the density of states in the momentum space, f{E) is as in (4) and F{E) 
is what was denoted by F{f) earlier, see (6). Performing one of the two integrations by 
using the delta function, we obtain 

W^^f = 271 J dE\ < t\H,nt\f > Wp{E)ff{E)F{E). (13) 

In order to evaluate the integral in (13), we note that 

f{E)F{E) = {f{E)f e(^-^-)/^, (14) 

where Ep is the Fermi energy. Notice that in the low-temperature limit we have 

f{E) = l/a E < Ef, 



and zero otherwise. Thus at low enough temperatures, the system for an arbitrary a (except 
very close to the Bosonic end) does exhibit a Fermi surface. We will be using this fact later 
in the calculation. Now for reasons that will become clear shortly, we wish to obtain an 
expression for the derivative df /dE. Using (4-5), we get 

^^---imf%---umr^, (15) 

Substituting {f{E)Y from (15) into (14), we get 

mF{E) = -|,Ii^e"=-')/-. (16) 

dtj uj[\ + uj) 

We now consider the low-temperature (T ^ Ep) limit of the process under consideration. 
In this limit, f\E) ^ 6{Ep — E)/a, and hence, df /dE ^ —6{Ep — E)/a. Substituting this 
in (16), we find 

f(E)F(E) « ^^^±^^S(E, - E). (17) 

Substituting (17) in (13) and performing the energy integration we get the Generalized 
Fermi Golden Rule: 



W,^f = 27r| < i\H,nt\f > ripiEp)) 



auj{EF){l + uj{EF))' 
In the special case when a = 1 (Fermions), we have, from (5), uj{Ef) = 1, and we get 

Wf_^f = 27r| < t\H,^t\f > \\p{EF)fT. 

We can therefore, in general, write 

'^f '-^faoo{EF){l + oo{EF)y 

We may now apply this result to specific cases. A straightforward application is to the 
calculation of a relaxation time of a nuclear spin in a metal. The relaxation time for arbitrary 
a is then given by*-^'^-' 

-{a) = —[hermions)- 



T T^ aij{EF){l+Uj{EF))' 
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Thus the change due to fractional statistics is simply given by a multiplicative factor which 
depends on a. In particular, at a = 1/2, this multiplicative factor can be explicitly calcu- 
lated, and we get 

— (a) = \^—{Fermions) . 

T T 

Notice that in deriving the Generalized Fermi Golden Rule we have made a number of 
simplifying assumptions. Real systems are likely to be more complicated. Nevertheless, the 
above derivation probably indicates the correct direction in which the transition rates move 
when fractional exclusion statistics particles are involved. 

In summary, we have generalized the Boltzmann-Nordheim equation and the Fermi 
Golden Rule for quasiparticles obeying the Haldane exclusion statistics. As two simpli- 
fied applications of these results, we have calculated the density response function and the 
relaxation time of a one-dimensional gas obeying the new statistics. Although the applica- 
tions presented here are simple and exploratory, we believe that they are the first attempts in 
this area to make a connection of theory to the consideration of observable effects. Moreover, 
this could be attempted only because the transport equation was nontrivially modified. 
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FIGURE CAPTION 

Fig. 1 : Dashed and solid lines represent integrals occurring in (10) and (12) respectively, 
as a function of eB/9/2, showing an enhancement of the density response function 
x{Q = 'ikp) at various temperatures, as one goes from Fermionic to semionic statistics. 
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